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Abstract 

The Dark Matter Time Projection Chamber (DMTPC) collaboration is develop- 
ing a low pressure gas TPC for detecting Weakly Interacting Massive Particle 
(WIMP)-nucleon interactions. Optical readout with CCD cameras allows for 
the detection of the daily modulation of the direction of the dark matter wind. 
In order to reach sensitivities required for WIMP detection, the detector needs 
to minimize backgrounds from electron recoils. This paper demonstrates that 
a simplified CCD analysis achieves 7.3 x 10~ 5 rejection of electron recoils while 
a charge analysis yields an electron rejection factor of 3.3 x 10~ 4 for events 
with 241 Am-equivalent ionization energy loss between 40 keV and 200 keV. A 
combined charge and CCD analysis yields a background-limited upper limit of 
1.1 x 10~ 5 (90% confidence level) for the rejection of 7 and electron events. Back- 
grounds from alpha decays from the field cage are eliminated by introducing a 
veto electrode that surrounds the sensitive region in the TPC. CCD-specific 
backgrounds are reduced more than two orders of magnitude when requiring a 
coincidence with the charge readout. 

Keywords: Dark matter, WIMP, TPC, CCD, Dark matter wind, Direct 
detection, Directional detection 



1. Introduction 

In recent years, dark matter direct detection experiments have obtained 
seemingly contradictory results, both supporting the existence of WIMP dark 
matter [lH3| and setting ever more stringent limits on its interaction cross section 
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with nucleons 0, HI ■ The tension between these experimental results highlights 
the need for a detection strategy that provides an unambiguous measurement 
capable of distinguishing between WIMP-nucleus scattering events and back- 
ground nuclear recoils. The strong expected directional signature of WIMP- 
induced recoils, due to the motion of the earth through the galactic dark matter 
halo, may provide such an unambiguous evidence of WIMP-nucleus scatter- 
ing [y, 0] • A number of experimental groups are developing detectors to search 
for this signal Q, and techniques have been developed to analyze data with 
directional information to extract the possible directional signal of WIMP dark 
matter (JQ. 



The DMTPC collaboration uses a low-pressure time projection chamber 
(TPC) to search for WIMP-nucleon elastic scattering from WIMPs in the local 
dark matter halo. The TPC is filled with CF4 gas to take advantage of the 
expected favorable WIMP- 19 F spin-dependent cross section [l5[. The WIMP 
interaction signature is a low-momentum nuclear recoil that leaves an ioniza- 
tion trail in the detector. Primary-ionization electrons from nuclear recoils are 



40 amplified 16| and the scintillation light from these avalanches is imaged by a 

41 charged-coupled device (CCD) camera. Using the shape of the ionization trail, 

42 DMTPC detectors are able to identify the direction of the nuclear recoil. The 

43 analysis of CCD tracks is described in (l7| . 

44 Electrons from j3 decays and processes such as Compton scattering from 7- 

45 and x-rays are typically important backgrounds in dark matter direct detection 

46 sources. These can be due to radioactive contaminants in the materials used to 

47 construct the detector, radioactive components of the target material, and ra- 

48 dioactive material in the laboratory environment. Such events are often rejected 

49 by fiducialization of the active volume. In experiments measuring the energy 

50 loss due to multiple processes, such as ionization, scintillation, and phonon ex- 

51 citation, these can further be identified by the relative fraction of energy loss 

52 due to each process 0, [B| . Other experiments attempt to create detectors that 

53 are wholly insensitive to electron recoils [3, . Many of these techniques are 

54 inapplicable to current directional detector designs, which typically measure 

55 only the ionization component of recoil energy loss. Directional detection ex- 

56 periments can instead use the low density of ionization from electron recoils to 

57 identify and reject these events 2(J 21 1. 

58 This work describes the charge readout systems of a prototype DMTPC 

59 detector and the corresponding analysis used to evaluate recoil properties such 

60 as position, energy, and geometry. It presents several studies taken with this 

61 detector in a surface laboratory at MIT. Section I2TT1 describes the charge readout 

62 systems of this detector and Section [3] describes event reconstruction. 

63 In DMTPC detectors, the CCD readout is used to reconstruct the ionization 

64 density projected onto a two-dimensional readout plane. The stopping power 

65 of electrons typically falls within the CCD read noise and is too low for these 
ee events to be reconstructed, while the stopping power of a particles and nuclear 
67 recoils is sufficient to accurately measure recoil properties. The charge signals, 
es which are sensitive to the total energy loss due to ionization, are able to measure 
69 and identify electronic recoils. Section |4~T1 describes a study using charge signals 
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70 to determine the ability of a DMTPC detector to reject electronic recoils and 

71 to enhance the ability of the detector to reject these events through pulse shape 

72 analysis. 

73 For DMTPC detectors, a number of backgrounds from CCD readout artifacts 

74 and interactions in the CCD bulk must also be eliminated in order to perform 

75 a low-background WIMP search. Section 14.21 describes these backgrounds and 

76 rejection strategies using the charge signals. It also discusses the use of a veto 

77 region to identify and reject alpha decays from the outside the sensitive region. 

78 In certain circumstances, a particles can be misinterpreted by the CCD as 

79 nuclear recoils. 



2. Detector Design 

The data in this paper were taken with a detector in a surface lab in 
Cambridge, Massachusetts using a 75 Torr CF4 gas target. Figure [T] shows 
a schematic of this detector. The field cage of the detector has a drift length of 
10 cm with a wire mesh cathode held at -1.2 kV. The drift cage is constructed 
of copper rings with a 27 cm inner diameter. The amplification region consists 
of a stainless steel woven mesh separated by nonconductive spacers from an 
anode plane made of copper-clad G10. The amplification mesh is coupled to an 
amplifier with a 30 f2 impedance to ground while the anode is biased at 680 V. 
The separation is 440 /im. 

The scintillation light from electron amplification goes through a viewport 
at the top of the vacuum vessel and into a Nikon f/1.2 lens with a focal length of 
55 mm. The lens directs the light into an Apogee Alta U6 CCD camera, which 
uses a Kodak KAF-1001 1024 x 1024 pixel CCD. The 24 ^m x 24 ^m pixels are 
binned on-chip into 4x4 blocks prior to digitization, resulting in a 256 bin x 256 
bin image. The CCD images a 16.7 cm x 16.7 cm square centered on the circular 
anode. The remainder of the anode is not imaged. Each 4x4 pixel bin reads out 
a 650 /im x 650 fim region of the anode. An americium-241 a source was was 
used to determine a CCD energy calibration of 13.1±0.1 analog-to-digital units 
(ADU) per keV of energy loss contributing to ionization. The energy of the a 
particles from the source is attenuated by a thin film, and the mean energy was 
measured to be 4.44 MeV using an Ortec ULTRA ion-implanted-silicon detector 



22[ . Calculations from SRIM [23J show that over 99% of the total energy loss of 
a particles at these energies contributes to ionization at this energy. Measured 
recoil energies in this paper are reported in a-equivalent energy units, denoted 
keV Q . Because the energy loss is dominated by ionization, the a-equivalent 
energy will be similar to the more standard electron-equivalent energy (in units 
of keV ee ), which is difficult to determine in a detector that is largely insensitive 
to electrons. 

2.1. Charge Readout Systems 

The anode plane has a diameter of 26.7 cm and is separated into two regions. 
The outermost 1 cm, called the veto region, is used to identify ionization events 
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Figure 1: A schematic of the detector: the drift field is created by a cathode mesh, field- 
shaping rings attached to a resistor chain and a ground mesh. The drift cage has a height of 
10 cm and a diameter of 27 cm. Primary ionization from a recoiling nucleus is drifted down 
to the ground mesh. The high-field amplification region is formed by the ground mesh and 
the anode plane. The grounded mesh is read out with a fast amplifier and the veto and anode 
are read out with charge-sensitive preamplifiers. The central anode region has a diameter of 
24.7 cm. Scintillation light from the amplification region is recorded with the CCD camera. 
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occurring near the field cage rings. The inner 24.7 cm diameter region, called 
the anode region, is used to measure the energies of recoil events. Charge 
signals from the anode region are amplified by a Cremat CR-113 charge sensitive 
preamplifier (CSP) Q, while the signals from the veto region are amplified by 
a Cremat CR-f 12 CSP. The gain of the CR-f 13 CSP is 1.5 mV/pC, and that 
of the CR-112 CSP is 13 mV/pC. Both have a nominal rise time of 20 ns when 
disconnected from the detector and a decay time of 50 [is. Recoil events have 
a typical rise time of approximately 1 /is due to the ion drift velocity across 
the amplification region, so the peak output voltage of the anode CSP gives a 
very accurate measurement of the total track ionization. Low-energy a particles 
from the 241 Am run described in Section I3T41 with 40 keV Q < E < 400 keV a , 
give an energy calibration of 0.251±0.003 mV per keV Q in the central anode 
channel. There is a linear relationship between the energies measured in charge 
and light in this energy range. 

The amplification mesh is read out through a Route2Electronics HS-AMP- 
CF preamplifier 25[. This preamplifier has a gain of 80 and a rise time of 
roughly 1 ns when not connected to the detector. A 30 ft resistor connects 
the mesh and preamplifier input to ground so that the output is proportional 
to the induced current from charged particles drifting in the amplification re- 
gion. Most ionization in the amplification region happens very near the anode 
so the electrons drift quickly over a very short distance while the ions drift more 
slowly over a longer distance, from the anode toward the amplification mesh. 
The induced current of the electrons from a single electron is a peak that decays 
within several nanoseconds while the ions create a broader but smaller shoul- 
der. Low-energy nuclear recoils create compact ionization trails, with a range 
along the drift direction (Az) of no more than a few millimeters. Because of 
the electron drift velocity of approximately 10 cm/ /is .26], all primary ioniza- 
tion electrons reach the amplification region within a period of several tens of 
nanoseconds. The electron signals of the avalanches add together to create a 
fast rising edge and a sharp peak in the current signal. The ions drifting in 
the amplification region then create a second, broader peak. This shape can be 
seen in Fig. 2(a) In tracks with large Az, such as most electronic recoils and 
minimum ionizing particles, the spread in time of primary ionization electrons 
entering the amplification region is long compared to the detector response of 
an electron avalanche, and the resulting pulse is characterized by a single broad 
peak from both the electrons and ions in the avalanche. Pulse shape analysis 
in the charge channels (see Fig. [2]) provides powerful discrimination between 
nuclear and electronic recoils and can strongly suppress such backgrounds in a 
rare-event search. 

All charge signals are digitized by AlazarTech ATS860 8-bit PCI digitizers 
using a sampling rate of 250 MHz. The digitizer bandwidth extends up to 65 
MHz. A total of 12288 samples (49.2 /xs) are saved with each trace, including 
4096 pre-trigger samples (16.4 /us). Charge events are triggered on the rising 
edge of either the mesh channel signal at 75 mV or the central anode signal at 
10 mV. This is sufficient to obtain a high expected efficiency for E >30 keV Q . 
During event readout the CCD is exposed for 1 s, while the digitizers collect 
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(a) Mesh (left), central anode (right, dark), and veto (right, light) signals from a 
75 keV a ol track 
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(b) 60 keVd e recoil, with signal from veto electrode 

Figure 2: Example smoothed traces of (a) an a track and (b) a vetoed electronic recoil. Note 
the two-peaked structure of the mesh pulse in (a), which is absent in (b). 



158 charge triggers. After each exposure the image and charge triggers (if any) are 

159 written to file for later processing and analysis. A camera shutter is not used so 

160 tracks occurring during the shifting and digitization of CCD pixels, a process 

161 taking approximately 200 ms, are seen in the CCD data but not in the charge 

162 data. 



163 3. Event Reconstruction 

164 3.1. CCD Tracks 

165 Recoil candidates are selected using the energy, range and several geometric 

166 moments computed from the CCD image. In order to maximize the recoil 

167 analysis efficiency, less stringent cuts were used than those described in [l7j |. 
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168 The CCD selection cuts are described in detail in Section |4~2"1 In Monte Carlo 

leg studies, the CCD track finding and cuts achieve over 70% efficiency for ionization 

170 yields greater than 40 keV Q and over 90% for ionization yields greater than 

171 50 keV 

172 3.2. Anode and Veto Signals 

173 The anode and veto signals are first smoothed using a Gaussian convolution 

174 with a — 80 ns. The typical rise time of a pulse in these channels is 1 fis, so 

175 the smoothing reduces the noise while having little effect on the pulse shape. 

176 The pre-trigger region is then used to determine the baseline voltage and noise 

177 RMS. The pulse shape is characterized by its peak voltage and time and the 

178 times on both the rising and falling edge where the pulse reaches 10%, 25%, 

179 50%, 75% and 90% of the baseline-subtracted peak height. 

180 3.3. Mesh Signals 

181 For the mesh channel, a Gaussian convolution with a = 6 ns is used to reduce 

182 the noise. After smoothing, the 25% to 75% rise time of the initial sharp edge 

183 of the electron signal in the data used in this paper is always greater than 15 ns. 

184 The broadening due to the smoothing here is no more than approximately 10%. 

185 Other features are less sensitive to the smoothing algorithm. 

186 As with the anode and veto channels, the baseline voltage and noise RMS are 

187 calculated from the pre-trigger samples. The initial electron peak of a nuclear 
las recoil candidate is identified as the first peak in a pulse with a height greater 

189 than 50% of the total peak height. The ion peak is defined as the largest peak in 

190 the pulse separated from the electron peak by more than 50 ns. The lowest point 

191 between the two peaks is also measured. The pulse shape is also characterized 

192 by similar rise and fall time variables as for the other channels. Rise times 

193 are calculated between the initial pulse baseline crossing and the electron peak 

194 and fall times between the ion peak and the baseline crossing at the end of the 

195 pulse. If only a single peak is found, rise and fall times are calculated from that 
we peak. Single-peaked events occur rarely because the reconstruction algorithm 

197 will typically identify a small noise fluctuation as an additional peak if no clear 

198 second peak is present. Fig. [3] shows an annotated example of a mesh pulse 

199 from an a particle. 

200 3.4- Nuclear Recoil Selection Criteria 

201 Selection criteria for nuclear recoils were determined by placing a 241 Am a 

202 source above the cathode mesh and outside the active volume of the detector. 

203 Most of the a energy is lost before crossing through the cathode into the active 

204 volume, so that a tracks of tens to a few hundred keV in energy are measured. 

205 Three sets of cuts are used to suppress background charge traces: (1) removal 

206 of triggers on noise and other pathological events; (2) removal of signals from 

207 tracks passing over the veto region; and (3) removal of events identified as 

208 electronic recoils (Table [T]). 
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Figure 3: Annotated mesh pulse of a 115 keV a a particle. A Gaussian convolution lias been 
used to reduce the noise prior to pulse shape characterization. The pulse start and end points 
are the nearest baseline crossings to the peak. Rise times are calculated from the electron 
peak and fall times are calculated from the ion peak. For anode and veto pulses, only a single 
peak is identified. 



Variable Description 

Va Peak anode channel voltage 

V\r Peak veto channel voltage 

25% to 90% rise time of the veto channel 

Tr 25% to 75% rise time of the mesh channel 

V e /Va Ratio of the mesh electron peak and the anode peak 

Vi/VA Ratio of the mesh ion peak and the anode channel peak 

Table 1: Description of reconstructed pulse shape variables used in analysis cuts. 

209 The electronic noise-reduction cuts remove events with anomalous baseline 

210 voltages or noise RMS voltages. Any charge triggers due to noise on one of the 

211 channels are removed using cuts on the pulse rise and fall times. The analysis 

212 also checks that the full mesh pulse is included in the saved waveform and that 

213 the mesh and anode pulses are correlated in time. Charge events that saturate 
2« the digitizer for one of the charge channels are removed. 

215 Events in the central anode region induce small pulses in the veto channel as 

216 well. Events actually passing through the veto region are rejected by requiring 

217 that Va > 4Vy and that < 400 ns. Pulses with longer rise times show 

218 the characteristic shape of electron avalanches occurring in the amplification 

219 region above the veto channel and are rejected with the latter cut. Finally, a 

220 small population of pulses showing pileup effects arc rejected by requiring time 

221 coincidence between the veto pulse and the mesh pulse. This population is very 

222 rare and typically only occurs when a radioactive source is placed inside the 
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Figure 4: Mesh signal rise time of both electronic recoils and a tracks. The dashed line 
represents the selection cut used in this analysis to remove electrons. 

223 detector. 

224 Nuclear recoil candidates are then identified using the shape of the mesh 

225 pulse. The rise time is due to a combination of electron longitudinal diffusion 

226 during drift, electronics response and recoil Az. Nuclear recoil candidates have 

227 < 22 ns for ionization yields of less than 125 keV Q (Fig. SJ. The compact 

228 tracks from nuclear recoils lead to generally larger but narrower peaks compared 

229 to electronic recoils of the same energy (Fig. [5]). In this analysis, nuclear recoil 

230 candidates have V e /VA > 4.5 and Vi/Va > 5.5. 

231 Once a set of candidate nuclear recoil light signals (using the cuts described in 

232 Section l4~2]) and a set of candidate nuclear recoil charge signals are identified, the 

233 event reconstruction attempts to match each track to its corresponding charge 

234 signal. To do this, all possible charge-light signal pairs are considered. The best 

235 match according to the relative charge-light energy calibration of KmodebriV] = 

236 (3.07 ± 0.04) + (0.01916 ± 0.00002)iVccD[AD[/] (Fig. [6]), where V^ode is the 

237 peak height of the anode signal and iVccD is the total number of ADU in the 

238 CCD track, is chosen. The match is accepted if the anode signal is less than 8.5 

239 mV (3.5cr) from the value estimated from the measured light signal. 

240 At the energies measured here, the trail of ionization left by a particles is 

241 much longer than that left by fluorine and carbon recoils of the same energy. The 

242 a particles in the calibration data enter the field cage at a mean polar angle of 

243 approximately 45° and drift across the full length of the drift cage (see the source 

244 placement in Fig. [1]), so the measured parts of the a tracks have more diffusion 

245 and Az comparable to or longer than nuclear recoils oriented exactly along the 

246 drift (z)direction. Because of this, the cuts set using the 241 Am data will be 

247 valid for carbon and fluorine recoils as well. The efficiency of the charge cuts 

248 and charge-light energy matching for tracks identified in the CCD analysis was 

249 measured using the a data to be greater than 90% for 40 keV Q < E < 200 keV Q . 



9 



> 

E 

o 

> 



137 Cs Y 
241 Ama 



n □ □ □ □ 
□ □ □ □□ 
u □□□ □□ 
□ □ □□ □□ 
° □ □ □□□□ 

□ □□□□□□ 
3 □ □ □ □□□□ 



□□□□□□□ O 

□ □ □ □ 

□□□□□□□□□□□□ □ □ 
■ □□□□□□□□□□□□ □ n 
o n a c 

innnanoDnan n □ □ n 




on/V a 



Figure 5: Ratio of mesh signal peaks to anode signal peak of electronic recoils and a tracks. 
The dashed line represents the selection cuts used to remove electrons. 
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Figure 7: 25%-75% rise time vs. i?2D of low energy a tracks. The tracks are collimated with a 



mean angle of roughly 45°, so Az increases with i?2D- The fit is to T(i?2D) = y^o m2 R\ 



where To = 17.19 ± 0.06 ns is the contribution to the rise time from diffusion and m = 
2.04 ± 0.03 ns/mm is the slope of the line at high -R2D- The error bars here represent the la 
spread from the mean value. A large spread is expected due to both straggling effects and the 
initial energy and angle distribution of the source. 

The a data was also used to evaluate the capability of using the rise time 
to measure the recoil Az. Because the mean polar angle of the a particles is 
roughly 45°, the two-dimensional range (i?2D) measured by the CCD is roughly 
proportional to Az. Due to straggling and imperfect collimation of the source, 
a range of a energies and incident angles is seen, so a wide range of Az values 
is expected for a given i?2D- However, the mean value can still be used to 
determine a calibration between the rise time and i?2D • Fig- El shows the result 
using the 25% to 75% rise time of the mesh signal. The proportionality constant 
between this rise time variable and i?2D is measured to be 2.04±0.03 ns/mm. 

4. Background Rejection 

4-1- Electronic recoils 

To determine the ability of a DMTPC detector to reject electronic recoils, 
a collimated 5 /xCi 137 Cs source was placed above the cathode mesh, outside 
the fiducial volume of the detector. The monoenergetic 660 keV 7 rays create 
a broad spectrum of electronic recoils across the entire energy range used for 
WIMP searches, while higher energy 512 keV and 1.2 MeV electrons from /? 
decays may also be detected. 

An average of 27 charge signals per 1 s exposure passed the noise and patho- 
logical event cuts during the Cs data run. Pile-up in the charge signals is 
negligible at this event rate. However, the signals of all recoils are accumu- 
lated into a CCD image, so that a region where several recoils overlap may be 
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Figure 8: Light: Energy spectrum measured by the central anode channel in the 13 'Cs run. 
Dark: Measured background spectrum scaled to match the live time of the 137 Cs run. The 
charge trigger efficiency decreases below 30 keV Q leading to the apparent peak in the spectrum. 



271 misidentified as a nuclear recoil, even if any single recoil would not be observable 

272 above the CCD read noise. 

273 A separate run with no sources inside the chamber was taken to measure 

274 the expected background spectrum. After removing events less than 3 seconds 

275 following a spark on the anode, a total of 51053 and 32568 one-second exposures 

276 were taken in the 137 Cs and background runs, respectively. The distribution of 

277 all charge events passing the initial set of noise and pathological-event cuts 

278 is used to determine the number of electronic recoils, as the trigger rate is 

279 dominated by electrons and minimum ionizing particles. The energy spectra 

280 of the two runs, after applying only the quality cuts identifying valid charge 

281 events, are shown in Figure [H 

282 In Monte Carlo studies, the CCD analysis used in this study has an efficiency 

283 of over 90% for nuclear recoils with energies between 40 keV Q and 200 keV a . In 

284 this energy range, a total of 679919 charge events were measured in the 137 Cs 

285 run with an expected background of 77310 ±350. In the CCD analysis 20 tracks 

286 in the background run and 63 in the 137 Cs run pass all cuts without considering 

287 any charge cuts. The higher event rate in the 137 Cs run provides evidence 

288 for possible light pileup events; many CCD tracks occur near the track-finding 

289 energy threshold. The number of charge traces passing each set of cuts is given 

290 in Tabled 

291 The camera only images part of the amplification plane while the charge 

292 channels read out the entire plane. A background-subtracted sum of images 

293 taken during the 137 Cs run is used to estimate the fraction of electronic recoil 

294 events occurring within the region viewed by the CCD. A fit of this sum to a 

295 two-dimensional Gaussian distribution shows that roughly 68% of all ionization 

296 from electronic recoils occurs in this region. This value is known to within 
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Table 2: Number of charge triggers with 40 keV a < E a node < 200 keV a that pass the specified 
cuts. The background subtraction includes the uncertainty of the background measurement. 

297 roughly 15%. It is expected that approximately the same percentage of recoils 

298 also occur in the region read out by the CCD. Using this number and subtracting 

299 the expected background, 409770±430 electronic recoils occurred in the fiducial 

300 volume out of a total of 602600±350 with 40 keV Q < E < 200 keV Q . Using 

301 these numbers, the CCD analysis alone achieves an electron rejection factor of 

302 7.3 x 10~ 5 in this energy range. The charge analysis includes all recoils occurring 

303 in the field cage and rejects electron recoils at the level of 3.3 x 10~ 4 . 

304 A combined analysis using both charge and CCD cuts as well as charge- 

305 light energy matching yields 5 events in the 137 Cs run and 3 events in the 

306 background run. Including the difference in total exposure time and the large 

307 uncertainties on both the background and combined background and signal 

308 distributions, a 90% confidence level upper limit on the electron recoil rejection 

309 factor of 1.1 x 10~ 5 is reached in the 40 keV a to 200 keVo, energy range of the 

310 recoil spectrum of 137 Cs 7 rays. From the detection efficiencies determined by 

311 Monte Carlo studies and the a source measurements, the detection efficiency 

312 for nuclear recoils using a combined CCD and charge analysis is determined 

313 to be over 63% for energies greater than 40 keV Q and over 80% for energies 

314 greater than 50 keV Q . A plot of the two-dimensional range as measured by the 

315 CCD against the recoil energy as measured by the anode channel shows (Fig. [9]) 

316 that the passing events appear near the predicted three-dimensional range for 

317 nuclear recoils from SRIM. The two-dimensional range is shorter than the full 

318 three-dimensional range, so nuclear recoils are expected in a broad band with 

319 shorter ranges than the SRIM prediction. Furthermore, the peak pixel values 

320 for the tracks are well above the threshold used in track finding, indicating that 

321 the events are likely to be nuclear recoil or a backgrounds rather than signal 

322 pileup from electrons. If these events were excluded as likely nuclear recoils, the 

323 result would be statistics-limited at 5.6 x 10~ 6 . 

324 Jf.,2. CCD Backgrounds 

325 Several classes of CCD background events appear in our detectors and have 

326 previously been eliminated through the use of CCD track selection cuts. While 

327 these cuts are generally able to identify and remove most such events, the charge 

328 readout analysis is also useful at eliminating these events. With the addition 

329 of the charge readout analysis, much less restrictive cuts can be used in the 

330 CCD analysis to maximize the efficiency for nuclear recoils. These background 

331 events include (1) residual bulk images, (2) intermittent hot pixels and (3) noise 
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Figure 9: 2-dimcnsional range vs energy of events in the 37 Cs and background runs that 
pass all cuts. Solid lines are SRIM predictions for carbon, fluorine and helium tracks. The 
measured range is two-dimensional, while the SRIM curves give three-dimensional ranges so 
the measured points will generally fall below the SRIM prediction. 



events. CCD backgrounds do not originate from ionization in the gas volume 
and so should have no associated charge signal. 

Residual bulk images (RBI's) generally occur in our detectors when a spark 
inside the amplification region causes a great deal of light to hit the CCD. A 
fraction of the longer wavelength light penetrates deeply enough into the silicon 
to generate photoelectrons in the depletion region. With little to no electric 
field in this region, the free charge carriers must diffuse via thermal motion into 
the pixel potential wells. This can take several minutes 27|. The strength of the 
RBI signal is proportional to the exposure time and can sometimes mimic low 
energy nuclear recoils. Residual bulk images occur in front-illuminated CCDs 
such as the Kodak KAF-1001 chip used in the Apogee Alta U6 camera but not 
in back-illuminated CCDs. 

Hot pixels and ionization from recoils occurring directly in the CCD chip are 
another source of background events found in the CCD event reconstruction. 
Hot pixels are a readout artifact where a single pixel is read out as containing 
an anomalously large amount of charge. Particles such as muons and 7 rays 
incident on the CCD can also leave ionization inside the CCD silicon that will 
appear in one or more pixels. Such tracks typically leave much more charge per 
pixel than light from avalanches inside the detector gas volume. 

While the cuts described in [TtJ are able to remove most of these CCD 
artifacts, they also can significantly reduce the reconstruction efficiency for low 
energy nuclear recoils, where nuclear recoils are more difficult to distinguish 
from the background noise of the CCD. Requiring coincidence between the light 
and charge signals provides additional rejection power of these types of events. 
The CCD analysis is less sensitive to recoils than the charge analysis due to the 
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Event type 


Before 


After 


Reduction [%] 


from CCD analysis 


charge cuts 


charge cuts 




RBI 


1332 


7 


99.5 ±0.2 


Hot pixel/CCD Si Recoil 


1246 


11 


99.1 ±0.3 


Edge Crossing (a) 


17 





iooi° 


Nuclear recoil/Out-of-time a 


20 


5 


75 ±11 


All Tracks 


2615 


23 


99.1 ±0.2 



Table 3: Number of CCD events of different types found before and after applying charge 
cuts, 25 keV a < -Eccd < 400 keV a . The RBI and hot pixel/CCD Si recoil events are false 
coincidences between low energy charge signals and CCD artifacts with low apparent light 
signals. 



357 favorable signal-to-noise ratio of the charge channels compared to the CCD. 

358 In addition to these CCD artifacts another CCD-related background is "out- 

359 of-time" events, which are recoils and a decays occurring during event readout. 

360 When the CCD is being read out, charge is shifted from one pixel to another as 

361 the charge from each pixel is digitized. A camera shutter is not used during data 

362 acquisition so any scintillation occurring during the 0.2-0.3 s readout time will 

363 appear shifted in the digitized image from its true position. The scintillation 

364 light from an a particle depositing only a small fraction of its energy in the 

365 imaged region can then be shifted toward the middle of the image, where it 

366 might resemble a nuclear recoil. These events are very difficult to identify with 
36? the CCD analysis. The charge channels do not collect data during this readout 

368 time, so a coincident charge signal will not appear for these shifted events. 

369 To evaluate the ability of the charge readout in removing these CCD back- 

370 ground events in the energy range 25 keV Q < Eqcd < 400 keVo,, a separate 

371 analysis was performed on the data from the background run described in Sec. 

372 14.11 A reduced set of CCD cuts identical to those used in Sec. 14.11 is used to 

373 define different classes of CCD artifacts and to determine the fraction removed 

374 by requiring charge-light coincidence. RBIs are defined as having at least two 

375 CCD tracks occurring within 10 pixels (1.6 mm) of one another within a single 

376 1000 event run. Hot pixels, ionization events in the CCD chip and noise events 

377 are identified by three cuts. They 1) have a maximum pixel value greater than 

378 500 ADU (38 keV), 2) include a pixel containing more than 30% of the total 

379 light of the reconstructed track, or 3) do not include enough pixels far enough 

380 above background to reconstruct a nonzero range. Tracks passing these cuts 

381 but touching the edge of the image are typically a decays and are removed as 

382 well. Any remaining tracks are tagged as nuclear recoil candidates or out-of-time 

383 partial a events. 

384 Requiring all CCD events to have a coincident charge signal removes more 

385 than 99% of the RBI, noise, CCD Si track, and hot pixel events before applying 

386 any of the CCD cuts. The few events of these types that did pass were due to 

387 false coincidences between CCD and charge signals with very low energies. All 

388 edge crossing events were removed, as would be expected since the CCD only 
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389 measures part of these events. Finally, 75 ±11% of the events tagged as possible 

390 nuclear recoils in the CCD-only analysis were also rejected, leaving 5 nuclear 

391 recoil candidates having both charge and light signals. This indicates that many 

392 of the events identified in the CCD analysis are likely either out-of-time nuclear 

393 recoils and a decays or events from the other classes that were not identified by 

394 the selection cuts (Table. [3]). 

395 5. Discussion and Conclusions 

396 This work has demonstrated the ability to reject electronic recoils by a factor 

397 of 1.1 x 10~ 5 at 90% C.L. level for electron recoils generated from a 137 Cs j- 

398 ray source with energies between 40 keV Q and 200 keV Q (or between 80 keV r 

399 and 300 keV r for fluorine recoils). Neither the charge nor the CCD analysis is 

400 completely efficient at removing electrons when considered independently in this 

401 analysis, but the combined result is significantly stronger. In an underground 

402 WIMP search, the CCD analysis is expected to be much more effective due to 

403 the greatly reduced recoil multiplicity per exposure compared to the 137 Cs run 

404 shown here. Stronger selection cuts will also further enhance the ability of the 

405 CCD analysis to reject electrons. The charge-light matching will also be much 

406 more effective as the chance of finding a false coincidence will be greatly reduced 

407 with the lower event rate in a source-free run. 

408 It is also expected that the CCD and charge analyses are most sensitive 

409 to different event topologies. The charge analysis is most effective at rejecting 

410 electronic recoils with large Az, while the CCD analysis is most efficient at 
4n rejecting electrons with large projected two-dimensional range, perpendicular to 

412 Az. With a rejection power in this simple analysis of 7.3 x 10 -5 in the CCD, 3.3 x 

413 10~ 4 in charge, and additional rejection power gained by requiring charge-light 

414 energy matching a rejection factor of order 10 -8 could be achieved if the CCD 

415 and charge analyses are in fact relatively uncorrelated. A much longer run would 

416 be required to attempt to confirm such a hypothesis. Even with this rejection 

417 power, the electron rejection analysis can be improved. Stronger cuts on recoils 

418 found in the CCD have already been used by the DMTPC collaboration in 

419 WIMP searches, although these also reduce the detector efficiency for signal 

420 events. Further refinements in the charge reconstruction and selection cuts are 

421 also expected in future analyses. 

422 The addition of the charge signal analysis strengthens the ability of DMTPC 

423 detectors to identify and remove CCD artifacts from the analysis by removing 

424 more than 99% of such events before applying CCD-based nuclear recoil selec- 

425 tion cuts. The charge analysis also allows for the identification and removal of 

426 events occurring during CCD readout, which can mimic the signal of low-energy 

427 nuclear recoils but cannot be easily identified with the CCD-only analyses used 

428 previously by the DMTPC collaboration. 

429 The ability to reject electronic recoils with high efficiency and without signifi- 

430 cantly reducing the detection efficiency of nuclear recoils suggests that electronic 

431 recoils are not expected to be a significant source of background events for the 

432 target energy range of DMTPC detectors for WIMP searches. 
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